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A free massive scalar field in inhomogeneous random media is investigated. The coefficients of 
the Klein-Gordon equation are taken to be random functions of the spatial coordinates. The case 
of an aiinealed-like disordered medium, modeled by centered stationary and Gaussian processes, is 
analyzed. After performing the averages over the random functions, we obtain the two-point causal 
Green's function of the model up to one-loop. The disordered scalar quantum field theory becomes 
qualitatively similar to a X(p^ self-interacting theory with a frequency-dependent coupling. 



I. INTRODUCTION 

Based on results obtained by Ford and collabora- 
tors [1, 2] and Hu and Shiokawa [3], recently three of 
us [4] proposed an analog model for fluctuations of light 
cones induced by quantum gravity effects. The model is 
based on two general features of waves propagating in 
random fluids. First, acoustic perturbations in a fluid 
define discontinuity surfaces that provide a causal struc- 
ture with sound cones. Second, propagation of acoustic 
excitations in a random medium are generally described 
by wave equations with a random speed of sound [5- 
8]. Taken together, these features lead to fluctuations of 
sound cones, analogous to the fluctuations of light cones. 
Phonons propagating in such a random fluid are then 
modeling photon propagation in a gravitational field with 
a fluctuating metric. 

In the present paper we depart from the quest for ana- 
log models of quantum gravity effects and study quan- 
tum flelds in a disordered environment in a more general 
context. Our study belongs to a wide program devoted 
to propagation of quantum matter fields in a classical 
background spacetime [9] but with metric fluctuations. 
Specifically, we consider a scalar quantum field described 
by a Klein-Gordon-like equation, in which the parame- 
ters of the equation, the mass and the coefficient of the 
second-order time derivative, become random functions 
of the spatial coordinates. The randomness of the pa- 
rameters is due to static noise sources that couple to the 
scalar field. In the limit of a vanishing mass, one re- 
covers the random wave equation considered in Ref. [4]. 
While not addressing the issue of the origin of the noise 
sources, we mention that one has in mind that they can 
be induced by a variety of phenomena, like metric fiuctu- 
ations due to quantum creation of gravitons in a squeezed 
coherent state in the presence of a black hole or interac- 
tions with background topological defects, among others. 
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Like in the case of quantum fields in the presence of an 
external heat bath, the noise sources break Lorentz sym- 
metry since they define a preferred reference frame. 

A scalar quantum field associated with acoustic waves 
in a disordered medium can define a situation where 
sound cones fluctuate randomly. In the study of such 
a situation, it is important to observe that systems with 
disorder can be divided into two wide groups, namely, 
systems with quenched or annealed disorder [10-12]. 
While in annealed systems the random field is in ther- 
mal equilibrium with the others degrees of freedom of 
the system, in quenched systems they are not. The dif- 
ferentiation between the two types of disorder is an im- 
portant issue in the studies of the infiuence of impurities 
on phase transition phenomena. Many authors have used 
field theory and the renormalization group to investigate 
systems with quenched disorders [13] , with particular in- 
terest in the role played by the disorder on critical expo- 
nents. In the case of quenched random fields, by means 
of the replica trick it is possible; to define a quenched gen- 
erating functional of connected n-point functions. With 
this generating functional in hand, the issue of the influ- 
ence of the disorder on phase transitions can be studied. 
In the present paper we consider annealed disorder. We 
consider weak noise fields and implement a perturbative 
expansion controlled by small parameters that character- 
ize the noise correlation functions. We obtain causal two- 
and four-point Green's functions of the scalar field. Aver- 
aging these Green's functions over the noise fluctuations, 
one obtains two- and four-point functions qualitatively 
similar to a self-interacting Xip"^ theory. As will be shown, 
we obtain a frequency-dependent coupling constant. 

Since Unruh's original paper [14], the possibility of 
simulating aspects of general relativity and quantum 
fields in curved space-time through analog models has 
been widely discussed in the literature [15, 16] . One inter- 
esting proposal is the generation of an acoustic metric in 
Bose- Einstein condensates and superfluids [17]. On the 
other hand, analog models with sonic black-hole could 
find a sort of generalization within the random fluid sce- 
nario. Since Bose-Einstein condensates are natural can- 
didates to produce an acoustic black-hole, one way to 
go beyond the semi-classical approximation is to inquire 
how these systems behave in the presence of disorder. 
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The calculations presented in the present paper are the 
first steps in the implementation of such a program. 

The organization of the paper is as follows: In Sec- 
tion II we discuss the perturbative approach in a free 
scalar field in the presence of annealed disorder described 
by two random functions. Section III contains our con- 
clusions. In the appendix some lengthy calculations are 
presented. 



II. PERTURBATION THEORY IN 
ANNEALED-LIKE DISORDERED MEDIA 



Let us consider a scalar field (p{t, r), defined in a (d+l) 
dimensional space-time, that satisfies the random Klein- 
Gordon equation: 



(1 + m) 



V2 + (l + 0m2 



ip{t,r)=0. (1) 



The quantities //, = /i(r) and ^ = ^(r) arc dimension- 
less random variables, functions of the spatial coordinates 
only, whose statistical properties will be defined shortly. 
Here, the random coefficient of the time derivative char- 
acterizes the possible discontinuity surfaces of the theory 
and provides a kind of random causal structure to it. 
It is well known that the principal part of a differential 
equation, i.e. the terms with the higher-order deriva- 
tives, determines entirely the loci of points in space-time 
where a solution may possess non-null discontinuities. 
This property is at the root of all analogue models of 
classical gravitation, namely, an effective causality can 
be obtained from the kinematical properties of a physi- 
cal system. The region of influence of excitations is given 
by an envelope which characterizes the maximum speed 
of propagation, that is, the light cone of the theory. From 
a physical point of view the region of influence of the the- 
ory may be obtained by an eikonal approximation, where 
the mass term may be disregarded. On the other hand, 
when the wave frequency is of the same magnitude of the 
mass term, excitations of the system propagate inside the 
characteristic cone. 

In the present paper we consider zero-mean random 
functions /u(r) and ^(r): 

(M(r))^ = 0, {ar)h = 0, (2) 

and, for simplicity, we suppose white-noise correlations: 

(M(rV(r')>^ = al6{r-r'), (3) 

(e(r)e(r'))e = a|<5(r-r'), (4) 

where {■ and (•••){ denote ensemble average of noise 
realizations and (7^ and cr^ characterize the strengths of 
the noises. We also suppose that the noises are Gaussian 
distributed (we use the notation fi{ri) = ^i): 

+ permutations, (5) 



and correlations of an odd number of noises arc zero. The 
reader may be aware, but it is worth remembering that 
these random functions fj, and ^ are statistically indepen- 
dent. 

Since Eq. (1) is linear in ip{t, r), it is useful to resort to 
Fourier transforms in order to find its solutions. There- 
fore, we define the Fourier transform of <f{t,r) on the 
time variable t: 

<p(t,r) = y^e--V('^,r), (6) 
and on the space variable r: 

^(•'^) = / (^^^'■""^(•''')- 

In addition, we define Fourier transforms of the station- 
ary noise functions: 



(8) 



and similarly for ^(r). Using these in Eq. (1), one obtains 
that the Fourier components of the field (^(a;,k) satisfy 
an algebraic equation that can be written as 



/ 



rfk' [io(k, k') + Li(k, k')] ^{cJ, k') = 0, (9) 



where Lq is a non-stochastic matrix with elements: 

Lo(k, k') = - k^ - mi^ d{k - k'), (10) 
and Li(k, k') is a stochastic matrix, with elements: 
Li(k, k') = ^ (^M(k - k') ^ - ^(k - k') mi^ .(11) 



In r-space, Lq and Li read: 



Loir) 



Li(r) = /i(r)-2 -^(r)m^. 



(12) 
(13) 



As in k-space, they act as integral convolution opera- 
tors. Note that while Li is diagonal, Lq is non-diagonal 
in r-space; yet, in k-space the situation is the opposite. 
In terms of these operators, the random Klein- Gordon 
equation can be written in matrix form: 



(Lo+ Li) (^(a;, •) = 0. 



(14) 



From this, one can define the full (operator valued) 
Green's function G: 



G={Lo + Liy\ 



(15) 
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Now, if one assumes that the noises are "weak" , a natu- 
ral perturbative expansion for G in the form of a Dyson 
formula can be defined: 

G = G(°) - G(°) Li G(°) + G(") Li G'") Li G'") + • • • 
= G(°)-G(°)EG(°), (16) 



with the self-energy S given by 



(17) 



where G^^-* = Lg ^ is the unperturbed (operator valued) 
Green's function which, in k-space, can be written as 



1 



- (k2 + ml) + ; 



(18) 



(Hereafter we take uq equal to unity.) A schematic repre- 
sentation of the expansion in Eq. (16) is shown in Fig. 1. 



FIG. 1: Perturbative expansion of G in terms of the disorder. 
The wavy lines represent generically the random functions fj, 
and 

We note that iG^°^ = Aq, where Aq is the noninter- 
acting Feynman propagator, which is the vacuum expec- 
tation value of the time ordered product of the quan- 
tum field operators at space-time points x = (t, x) and 
x' = (t',x'): 



Ao(x,x') = (0|r[(p(t,x)^(i',x')]|0). 



(19) 



Let us write down explicitly a few terms of the pertur- 
bative series in r-space, up to second order in the random 
fields: 



G{uj, r, r') = G(") (a;, r, r') - J dn G^°\oj, r, n) [^(ri) cj' - ^ri) ml] G^°\oj, n, r') 

+ JdnJ dr2G^°\oj,r,r2)[^i{v2)oj^-ar2)ml] G(°)(c^, r2, ri) 
X [f,{v,)oj'-ari)ml] G(°)(c.,ri,r') -(-••• . 



A pictorial representation in r-space of this perturbative 
expression in terms of diagrams can be done. They cor- 
respond to diagrams of multiple-scattering of the Klein- 
Gordon field on random inhomogeneity scatterers located 
at positions ri, r2, • • • . In k-space, similar diagrams cor- 
respond to multiple interactions of Fourier components 
of the Klein-Gordon field and of the random inhomo- 
geneities. 

Let us now perform the averages over the random pro- 
cesses that appear in the definition of the propagator G 
given in Eq. (20). Note that, because of the Gaussian na- 
ture of noise correlations, terms with an odd-number of 
noise fields do not contribute to the two-point function, 
and the first nonzero correction to the two-point function 
comes from the averages over the third term in Eq. (20). 
Therefore, up to second order in the noise fields, we get 

GW(a;,r,r') ^ ( G(c., r, r'))^^ 

= G(°)(^,r,r') + GW(w,r,r'), (21) 

with the one-loop correction G*^^-' (w, r, r'), represented 
pictorially in Fig. 2 and given in terms of its Fourier 
transform G'^-'(a;,k) as 

GW(^,r,r') = J ^G(i)(u;,k)e^'^-('--'). (22) 

This quantity can be written in terms of a self-energy 



(20) 
I 

S(a;,k) - defined consistently with Eq. (16) - as: 

G(i)(w,k) = -G(°)(cj,k)I](cj,k)G(°)(w,k), (23) 

with 

E(c., k) = - {ay + aim',) a{co), (24) 

where 

(l^^^°^('"'^)- (25) 

Note that at one-loop order, the self-energy is actually in- 
dependent of k. From Eq. (16), one has G~^ = G*^°^~-^-|-E 
and, therefore, in the k-representation 

[G(i)(w,k)]-i = [G(o)(c.,k)]-i + I](c.,k) 

= uj^ - m^, (26) 

where 

— mQ — I](aj,k) 

= ml + {ay + almt) a{uj). (27) 

From this result, it is clear that the effect of ran- 
domness is to turn a free, conventional scalar quantum 
field theory into a self-interacting theory, with a self- 
interaction qualitatively similar to Xip"^. The induced Xip'^ 
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FIG. 2: The one-loop correction to the two-point causal func- 
tion due to the random inhomogeneities. 



theory is a consequence of the Gaussian nature of the 
functions ^(r ) and ^(r); a more general polynomial self- 
interacting model is obtained if non-Gaussian noises are 
used. The integral is divergent because of the white-noise 
nature of the noise; a colored noise correlation function 
would lead to a finite integral. 

As mentioned earlier in the text, the induced coupling 
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natural to expect that such function can bring us some 
information on the nature of the induced coupling. 

Since the random field equation is linear in the field 
variable, the corresponding action of the system is 
quadratic in the field, resulting in a Gaussian generating 
functional. Then, the quantum rt-point Green's func- 
tions of the system are products of two-point Green's 
functions. In particular the four-point Green's function 
will be of the form: 



G4{X1,X2,X3,X4) 



{G{xi,X2)G{x3,X4))f_,^ 
{G{xi,Xs)G{x2,X4))^i: 

{G{xi,x4)G{x2,X3))f_,^. (33) 



As before, x = (i,x). Fig. 3 presents a graphical repre- 
sentation of this expression. 



generated is frequency-dependent, A(w) ^ af^uj' 
Also, for c? = 3 one can isolate the finite part from the 
integral making use of the identity 



|k|2 _ (a;2 -to2) 



1 



uj — rrin 



|k|2 - (a;2 - mg) 



(28) 



so that the finite part of the self-energy, (w, k) , is given 
by 



Ey(c.,k) = (aX + a|m4) / 



d|k| uj"^ — rriQ 



27r2 w2 _ |k|2 _ 



FIG. 3: The four-point Green function up to one-loop, after 
noise averaging. 'IPR graphs' means one-particle reducible 
graphs. 

In Appendix A, Eq. (33) is computed up to one 
loop; the result can be written as G4(a;i, a;2, X3, X4) = 



G^^\xi,X2,X3,xa) + G^l' {xi,X2,X3,Xi), where G^^' and 



(0) 



^ (fj^w^ -I- a^m*) |a;^ m^l^/^yl (29) G4^'' are the tree-level and one-loop contributions, respec- 

noise averaffinff leading to the tree-level four- 



wit h 



A 



rriQ > w 



0- 



(30) 



tively. The noise averaging leading to the tree-level four- 
point function Gf'' is illustrated in Fig. 4 and the one- 
loop corrections gI^^ are pictorially represented in Figs. 5 
and 6. 



The conclusion from this calculation is that the ran- 
dom fluctuations induce a self-energy with a width, which 
gives a finite life-time for the excitations. 

Since randomness, as mentioned previously, led to an 
induced coupling, one may ask on the form of the per- 
turbative corrections to the four-point function G4. The 
reason for this particular interest is the following. As 
well known, given the four-point function, one can define 
the one-particle-irreducible (IPI) four-point connected 
Green's function without the external legs r4: 




r4(fci,fc2,fc3,fc4) = G4(fci,A:2,fc3,fc4) JjFfe,, (31) 

where 



i=l 



(32) 



In a covariant quantum field theory, one may compute 
the renormalized coupling constant from the proper ver- 
tex function, under certain conditions. So, even though 
we are not dealing with a covariant field theory, it is 



FIG. 4: The disorder-induced, tree-level four-point Green's 
function. 



From Eq. (31) and the results derived in the Appendix, 
one can express the tree-level proper vertex function as: 

Tf \ki,k2,k3,k4) = {2tt) (5 (^3 + ^4) (ct^W2^4 + (t|too) 
+ 6{u}2+ W4) {(^l^l^l + o-fmo) 
+ S{u:2+i^3){(ylujlujl + alml)\.{M) 

As mentioned earlier, the proper understanding of the 
nature of the induced interacting theory requires the 
study of the four-point Green's function. We will not 
dwell into such a study here. Besides being involved (see 
the Appendix for the one-loop contribution), such a study 
is out of the scope of the present paper and therefore we 
reserve it for a future publication. 
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FIG. 5: One-loop corrections of first kind to the four-point 
Green's function. 



The most striking feature that comes out from these 
calculations is that the induce coupling is qualitatively 
similar to a theory, with a frequency-dependent cou- 
pling. The important meaning of this is that if one starts 
with an interacting model, the induced coupling due to 
the random scatterings in the medium changes the value 
of the renormalized coupling constant of the model. As a 
result, a relevant question to be answered is to know what 
would be the sign of this effective coupling. For instance, 
if it turns out to be negative, it is clear that we would 
have a decrease in the value of the renormalized cou- 
pling constant. Furthermore, as well known, a negative 
coupling constant can be a source of instability in field 
theory, which means that one should look for alternative 
solutions to circumvent the above mentioned problem. 
One possibility is to discuss the ground state properties 
of the system in a theory with metastable vacua. This 
can be obtained in non-simply connected manifolds. For 
simplicity one might assume periodic boundary condi- 
tions in all d + 1 directions, with compactified lengths 
Li, L2,..,Lii^i. At this point we remark that in theo- 
ries defined on a non-simply connected Euclidean space, 
one of the compactified length is related to the temper- 
ature, remembering the Kubo-Martin-Schwinger (KMS) 
condition [18]. Therefore in a Euclidean theory describ- 
ing bosons, we have to impose periodic boundary con- 
dition in one compactified direction. In the other d di- 
rections, we are free to choose any boundary conditions. 
However, imposing periodic boundary conditions in all 
compactified dimensions enables us to maintain trans- 
lational invariance [19]. Non-translational invariant sys- 
tems were studied in Refs. [20]. We take this opportunity 
to call attention for the results obtained by Arias and co- 
workers [21], where the thermodynamic of the massless 
self-interacting scalar field model with negative coupling 
constant was analyzed. We believe that such results from 
this reference can be useful to proceed with the investi- 
gations of the consequences of impurities in a relativistic 
model. 




FIG. 6: One-loop corrections of second kind to the four-point 
Green's function. 



III. CONCLUSIONS 

Recently it was proposed that an analog model for 
quantum gravity effects can be implemented in condensed 
matter physics when the fluid is a random medium [4] . In 
order to study the wave propagation of the elementary 
excitations in the fluid when the acoustic wave propa- 
gates in random medium, it was considered that there 
are no random fluctuation in the density, but only in the 
reciprocal of the bulk modulus. In this analog model 
of quantum gravity with a colored noise a free quantum 
field theory describing phonons becomes a self-interacting 
model. 

In this paper we studied a massive scalar field theory in 
a disordered medium. We determined the perturbation 
theory in annealed-like disordered medium, with two dis- 
tinct random functions. After performing the random av- 
erages over the noise function a free scalar quantum field 
theory became the A(/?^ self-interacting model. The one- 
loop two and four-point functions were presented. We 
obtained that the induced coupling constant is frequency 
dependent. Note that although we employed Gaussian 
functions, it is not difficult to extend the method to non- 
Gaussian functions by means of a cluster expansion [22] . 

There exist several possibilities to use the above results 
in analog models. For example, a random fluid with a 
supersonic acoustic flow can be an analog model that 
allows one to discuss the effect of the fluctuation of the 
geometry in the Hawking radiation. Since Bose-Einstein 
condensate is a natural candidate to produce acoustic 
black- holes, in order to introduce randomness, i.e., the 
effects of the metric fluctuations in the acoustic black- 
hole, we have to investigate how superfluid and Bose- 
Einstein condensate behave in the presence of disorder. 
Sound waves in Bose-Einstein condensate, described by 
random wave equation must reproduce physics beyond 
the semi-classical approximation in an analog model. 

An interesting course of action would be the search of a 
possible localization of the Hawking thermal flux in such 
analog model, since it is well known that the effect of 
the impurities is the localization of classical waves and 
elementary excitations [23]. In conclusion, the study of 
quantum fields in disorder medium in the analog mod- 
els scenario introduces new experimental and theoretical 
challenges. The study of the relativistic Bose-Einstein 
condensation in the presence of a random potential [24] , 
and also the consequences of introducing metric fluctua- 
tions in the analog model proposed in Ref. [25] is under 
investigation by the authors. 
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Appendix A: Detailed derivation of the causal two- 
and four-point functions 

We start from the perturbative expansion for the two- 
point function G, Eq. (16), which we rewrite as: 

G = G(") - G(") Li G(") + G("' Li G^"' Li G^") + • • • 
with Q^'^^^ given by 

In terms of time and space variables, x — (t,x), this 
corresponds to 



Gix,x') = I dziG("^(a;-zi) 
and 



(A3) 



g(")(zi,.T')=(-i)"n^i(^j-) y rf^j-+iG(°n^.-,^.-+i), 

(A4) 



where ^1(2:) is given by 



Li(x) = Li(i,x) = -Ai(x)^ - C(x)m2. (A5) 

In terms of the frequency and wave-number coordinates 
in Fourier space, (w,k), Li is given by Eq. (11). In 
Eq. (A4), it is to be understood that z^+i = x' and 
that there is no integration in Zn+i- In Fig. 7 we present 
a graphical representation of a generic term in equation 
(A4). 



X=Zn+i Zn "' Z2 Zi X=Zo 



FIG. 7: Graphical representation of a generic term in Eq. (A4) 

For completeness, we give here the expression for 
Fourier transform of the free Green's function G'^^^ {x — 
x'): 



G(°)(a;,x') = 



dk 



(27r)" 



_g-«fe(x-:.')(^(0)(;.)^ (A6) 



with 



G(°)(fc) = G(°)(c^,k)=r^i 



1 



(A7) 



From the properties of the noise functions, given in 



Eqs. (2)-(4), one has that: 

(Li(x))^4=0, 



and 

(Li(x)Li(x'))k = [cjI^ ^+^lrnt^ ^(x~x'), 

(A9) 

From this, it is clear that: 

{g^^-+^\x,x')),^^Q, (AlO) 

and, therefore, the first correction to the two-point causal 
Green's function is given by: 

G^^\x,x') = Jdz,G^'>\x,z,){g^'\z,,x')),,^. (All) 

Now, for n — 2, Eq. (A4) gives: 



{g^^^izi,x'))^^ = I dZ2{Lliz,)Ll{z2)) 



X G(°)(zi,z2)G(")(z2,a;'). (A12) 

Muhiplying this by G'-"\x - zi ) and integrating over zi, 
one easily obtains the result: 



G'^^\x,x') 



dk 



_^-^k(x-x')GW(^k), (A13) 



with 



(A8) 



= -G(°)(^,k)I](w,k)G('')(cj,k), (A14) 

where T,{uj,'k) is the self-energy given in Eq. (24). 

Next, we derive the four-point function. Since the full 
action of the model is quadratic in the field, the quantum 
four-point function can be factored into the product of 
two two-point functions. The average over the noises of 
this function can then be written as 

G4:{xi,X2,X3,X4) = {G{xi,X2)G{x3,X4))f^^ 
+ {G{xi,Xz)G{x2,Xa))^^ 

+ (G(xi,a;4)G(a;2,X3))^e.(Al5) 

We are interested only in the one-particle irreducible 
(IPI) parts of this. The lowest-order IPI contribution 
is the tree-level four-point function - see Fig. 4. For first 
term in Eq. (A15) one has: 

(G(a;i,X2)G(x3,X4))^°^ = j dz^ J dz2G^"Hx,, z,) 

X G(o) (X3, Z2)( a^^' {zi,X2)g^^\z2,X4) (A16) 
For n ~ 1, Eq. (A4) gives: 

(gW(zi,a;2)e<'H^2,a;4))^5 = (Li(zi)Li(z2))^e 

xG(°)(^i,^2)G(o)(^2,a;4)(A17) 

Using the result of Eq. (A9) and performing the appro- 
priate Fourier transforms, one obtains: 
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dki f dki 



[att)"-^^ J (27r)"+-'- 

X ((T^cj^^l + alm.l) {2nY+H{ki + ^2 + ^:3 + A:4)(27r)(S(a;3 + oja). (A18) 

To arrive at the last result, we have made the changes k2 — > —k2 and k^ — >■ — ^4, and used the identity (5(ki + k2 + 
ks + k4)(5(wi + (jj2)5{ujz + W4) = 5{ki + ^2 + ^3 + k4)5{(jJz + ^4)- The other tree-level contributions in Eq. (A15) can 
be obtained from the result in Eq. (A18) by appropriate changes of integration variables. Specifically, the final result 
for the tree-level four-point function, Gf\xi,X2,x^,X4), can be written as: 

Gr(.„.2,.3,.4) = / J J^{2.r'Sik,+k2 + ks + k,) 

^ g-i(fclXl-Ffe2X2-Ffe3X3-Ffc4^4)C!W(fc^^^2,fc3,fc4), (A19) 

with 

Gf (fci,fc2,fc3,A;4) = G(0)(A;i)G(°)(^2)G(0Hfc3)G(°)(A;4)rf (A;i,fc2,fc3,fc4), (A20) 

where 

Tf\ki, fc2, fcs, k4) = {2tt) 6{uj3 + UJ4) {a^^ujlujl + CT|mo) 

+ 5{u)2 + W4) ("■^'^3^1 + (T|mo) + 5{(j02 + W3) {(^ft^l^l + i^l^^o) • (■'^^l) 

Next, let us consider the one-loop corrections to the four-point function. Let us focus on contractions which are 
depicted in Fig. 5. For the first term in Eq. (A15) we have: 

(G(a;i,X2)G(a;3,a;4))^^j"^^ = j dzi Xi, Zi)G^°^ {Xs, Z3)( G''^\zi,X2)G'''^\z3, X4) (A22) 

From Eq. (A4) , one has: 

{G^^\zuX2)g^^\z3,X4))^,i = j dz2 j dZ4{Li{zi)Li{z2)Li{z3)Li{z4))f,i 

X G^°\ZI,Z2)G^°\Z2,X2)G^°\Z3,Z4)G^'>\Z4,X4). (A23) 

Using now the fact that the noise correlations are Gaussian : 

(Li(zi)ii (2:2)^1 (z3)Li(z4))mC = {Li{zi)Li{z3))^i {Li{z2)Li{z4)} + {Li{zi)Li{z4)) (Li (22)^1 (z3))a.5(A24) 

The contractions {Li{zi)Li{z2))ij,^ (Li(^3)Li(z4))p,^ do not appear above since they do not contribute to one-loop 
corrections to the four-point function. Collecting the above results and using Fourier transforms, we have the first 
contribution to Eq. (A22): 

{G{x„X2)G{x3,X4))%-"'^ = J J0^---J ^^^e-^('=^--'==-^+'=3X3-fc4.4)G(o)(fc^)G(o)(fc2)G(°)(A:3)GW(fc4) 

X G(°)(fc5)G(°)(fc6) (->2'-4' + <ylmt) {al^l4. + 4^^) 

X j dZl--- j dZi e*('=l-'==)^l+»('==-'=2)^2+*('=3-fe6)^3-Ki(fc6-fe4)z4 ^(21 - Z3)^(Z2 - Z4) 

f '^^6_g_i(fe,,i+fe.x.+fc3X3+fc4X4)cr(0)(^.^)(;.(0)(^^)(^(0)(^^)g(0)(^^) 



{2ttY+^ j (27r)'^+i 
X G(°)(fc5)G(°)(fc6) {<yl^W4 + 40 K^We + 40 (27r)^+M(fci + fc2 + fc3 + ki) 
X {2'nY+H{k2 + k4 + kr, + kQ){2'K)6{L03+uj4){2T^)6{oJ2+oj^)- (A25) 
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To reach the last result we used the identities: 

(5(-ki - ks + ks + k6)5(k2 + k4 - ks - kg) = <5(-ki + ks - kg + k4)^(k2 + k4 - kg - kg), (A26) 

= 5{loi — UJ2 + 00s — u)4)S{u;s — u)4)S{—co2 — C04 + 1^5 + i^6)S{i^5 — 1^2), (A27) 

and, finally 

S{u!i -UJ2+UJ3- uJ4)S{uj3 - UJ4)6{-UJ2 - CO4 + LO5 + uj6)5{lj5 - uj2)S{-'ki + k2 - ks + k4)(5(k2 + k4 - k5 - ke) 
= S{-ki +k2 - k3 + k4)5{k2 + ^4 - fcs - kQ)6{(jj3 — a;4)5(a;5 - ^2) (A28) 

Then, we have made the changes k2 —k2 and ki —ki. 
As for the second contribution to equation (A22), we have: 

{G{x„X2)G{x3,X4))^;f"'^ = I (^^•••/ ^^^e-^('^--'=-^+'=3-3-fe.-4)G(o)(fc^)G(o)(fc2)G(°)(fc3)G(°)(A;4) 

j dZi--- j d^4e*('=i-*5)^i+^('=^-*=)^=+*(*3-*6)^3+i(fc6-fc4)^4^(2.1 - Z3)5(Z2 - Z4) 



X 



(27r)''+i J i27rY+i 
X G(°)(fc5)G(°)(fc6) {cjlioiui + afm^) (a^w^wl + a|m4) (2^)'^+M(fci + ^2 + ^3 + ^4) 
X (27r)'^+i5(fc2 + fc4 + A:5 + A:6)(27r)<5(w3+W4)(27r)(5(w2+W5). (A29) 

The last step in the above equation was achieved through the identities: 

(5(-ki + k4 + ks - k6)<5(k2 - ks - ks + kg) = <5(-ki + k2 - ks + k4)<5(k2 - ks - ks + kg), (A30) 

S{lOi - LJ5)S{U5 - LJ2)S{l03 - UJe)5{LUij - LU4) = S{uJi - UJ2)S{UJ3 - UJ4)5{lJ5 - iJ2)S{LJ3 - We) 

= -a;2 +ws - W4)5(a;s - W4)5(-W2 +ws +ws - W6)(5(ws - W2), (A31) 

and, finally 

S{lji -UJ2+0J3- 0J4)6{uj3 - UJ4)S{-UJ2 + UJ3 + LJ5 - uj6)S{lj5 - uj2)S{-'ki + k2 - ks + k4)(5(k2 - ks - ks + ke) 
= S{-ki +k2- k3 + k4)5{k2 - fcs - fcs + kQ)5{ui3 - uj4)5{uj5 - UJ2) (A32) 

Similarly as before, we have also made the changes k2 — )• —k2, ^4 — )• —k4 and fee — fee- 

The second kind of contractions that also lead to one-loop corrections to the four-point function is shown in Fig. 
6. For the first term in Eq. (A15), we have: 

(G(a;i,a;2)G(a;s,a;4))^j~^^^ = J dzi J dZ2G^°\xuZl)G^''\x3,Z2){g^^\zi,X2)g^^\z2,X4))^,S. (A33) 

From Eq. (A4), one has: 

{g^^\zi,X2)g^^\z2,X4))^,^ = j dz2 j dz4{Li{zi)Li{z2)Li{z3)Li{z4))^i 

X G(°)(zi,X2)G(°)(z2,Z3)GW(z3,^4)G(°)(^4,X4). (A34) 

Using again the fact that the noise correlations are Gaussian and since we are looking for one-loop corrections, we 
get: 

(Li(2i)Li(z2)Li(z3)Li(z4)),.4 = (Li(zi)Li(z3))^j (Li(z2)ii(^4));,5 (A35) 
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The other contractions of the decomposition above do not appear since they do not contribute to one-loop corrections 
to the four-point function. CoUecting the above results and using Fourier transforms, we have: 



{G{xi,X2)G{Xi,XA))\^ 



(1-//) 



dki 



dkA 



dk 



dk' 



)GW(fci)GW(fc2) 



(27r)'^+i (27r)^+i J {2t^Y+^ {2ttY+^ 

dZl--- j dZ4 e*('=l-fe)^l+'('=3-fc)^2+i(fc-fc')^3+i(fe'-fc4)^4 - Z3)5(Z2 - Z4) 

dki dki f dk dk' 



-I 



—i{kixi+k2X2+k3X3+k4,X4,)r'(0) (I. \r<{0)(i. ^ 

{2Try+^ {2Try+^ J {2ttY+^ {2ttY+^ ^ 

X S{ki + k2 + k3 + ki) {2'Kf+'^S{ki + k2 + k + fc')(27r)<5(a;i + a;2)(27r)^(a; + lj'). (A36) 

In the last equation we used similar identities of delta functions which were employed in calculation of the expressions 

(A25) and (A29). Just as the tree- level case, all the one- loop contributions in equation (A15) can be evaluated from 
the results in equations (A25), (A29) and (A36) by adequate interchanges of variables. Thus, the final result for the 
one-loop corrections for the four-point function, gI^\xi,X2,x^,X4), reads: 

G^^\xuX2,x„x,) = I j^-'-J j^{^-r'S{k, + k2 + ks + k,) 



with 



where 



X ^-i{k^x^+k2X2+ksxs+k^x^) G^^ [kxM^kzM), 



Gi'\ki,k2, fc3, ki) = G^°\ki)G^'>\k2)G'-°\k3)G^°\ki) T^PikiMM, ki), 



(A37) 



(A38) 



r^^\ki,k2,k3,ki) = (27r) 



d+3 



dk 



dk' 



YGW(fc)C?(°)(fc') 



(27r)'i+i J {2TTy 
5{k2 + ki^- k + k')5{uj3 + uj4)5{uj + uj2 



+ 6{k2 + k3 + k + k' 

+ S{k3 + k4 + k + k' 

+ 5{k2 + k-i + k + k' 

+ 6{k3 + k4, + k + k' 

+ 6{k2 + ki + k + k' 

+ 6{ki+k2 + k + k' 

+ 6{ki+k3 + k + k' 

+ 5{ki +k4 + k + k' 

+ 6{k2 + k3 + k + k' 

+ 6{k2 + ki + k + k' 

+ 5{k3 + ki + k + k' 



)5(u]3 + U!4)5{u] + LO2) 

)5{uj2 + oj4)5{uj + ^3) 

)5{UJ2 + i04)5{u) + UJ3) 
)6{l02 +U13)5{UJ -I-W4) 
)5{(jJ2 + W3)^(w -I- W4) 
)5{lOi + LJ2)^{i^ + w') 
)5{lOi + ijJ3)5{LO + LO') 
)d{iOi + UJ4,)d{iO + Lo') 
)S{U)2 +UJ3)S{U) +U)') 
)S{C02 + L04,)S{cO + Co') 
)d{i03 + L04)d{L0 + LO') 



( 2 2 2, 2 4\/'22/2i 24\ 



2 2 2 2 4\ / 2 /2 2 

2 2 2 I 2 4\ / 2 2 /2 I 2 4\ 

^^^^^3^4 + o'c'^o) ^ + o-c»™o) 

2 2 2 I 2 4\ Z 2 /2 2 I 2 4\ 

222, 2 4\ / 2 2 /2 , 2 4\ 

2 2 2 I 2 4\ / 2 /2 2 I 2 4\ 

2 2 /2, 2 4\ / 2 2 2 , 2 4\ 

2 2 /2, 2 4\ / 2 2 2 , 2 4\ 

cr^W3W + (T^mo) (ct^w W4 -I- a^mo) 

2 4\ / 2 2 2 , 2 4\ 

2 4\ Z 2 2 2 , 2 4\ 

2 4\ / 2 2 2 , 2 4\ 

2 2 /2, 2 4\ / 2 2 2 , 2 4\ 



,2™4\ 



2 4\ 



2 2/2 
2 2/2 
2 2/2 



(A39) 



10 



[1] L. H. Ford, Phys. Rev. D51, 1692 (1995). 

[2] L. H. Ford and N. F. Svaitcr, Phys. Rev. D56, 2226 
(1997); L. H. Ford and N. F. Svaiter, Phys. Rev. D54, 
2640 (1996); H. Yu and L. H. Ford, Phys. Rev. D60, 
084023 (1999); R. T. Thompson and L. H. Ford, Phys. 
Rev. D78, 024014 (2008); R. T. Thompson and L. H. 
Ford, Class. Quant. Grav. 25, 154006 (2008); H. Yu, 
N. F. Svaiter and L. H. Ford, Phys. Rev. D80, 124019 
(2009). 

[3] B. L. Hu and K. Shiokawa, Phys. Rev. D57, 3474 (1998). 

[4] G. Krein, G. Menezes and N. F. Svaiter, Phys. Rev. Lett. 
105, 131301 (2010). 

[5] A. Ishimaru, Wave Propagation and Scattering in Ran- 
dom Media (Academic, Now York, 1978). 

[6] H. Gzyl, Diffusion and Waves (Kluwer Academic Pub- 
lishers, Netherlands, 2002). 

[7] V. I Klyatskin Dynamics of Stochastic Systems (Elsevier 
B. v., Amsterdam, 2005). 

[8] J. Pierre Fouquc, J. Garnicr, G. Papanicolaou and K. 
Solna, Wave Propagation and Time Reversal in Ran- 
domly Layered Media (Spinger Science and Bussiness Me- 
dia, New York, 2007). 

[9] N. D. Birrell and P. C. Davis, Quantum Fields in Curved 
Space (Cambridge University Press, New York, 1982). 
[10] D. Chowdhury, Spin Glasses and Other Frustrated Sys- 
tems (World Scientific, Singapore, 1986). 
[11] M. Mezard, G. Parisi, and M.A. Virasoro, Spin Glass 

Theory and Beyond (World Scientific, Singapore, 1987) 
[12] G. Parisi, Lett. Math. Phys. 88, 255 (2009), S. W. Ed- 
wards and P. W. Anderson, J. Phys. r4, 965 (1975), 
D. Sherrington and S. Kirkpatrick, Phys. Rev. Lett. 35, 
1972 (1975). 

[13] T. C. Lubcnsky, Phys. Rev Bll, 3573 (1975), G. Grin- 
stein and A. Luther, Phys. Rev B13, 1329 (1976), J. A. 
Hertz, Phys. Rev. B18, 4875 (1978), S. K. Ma and J. 
Rudnick, Phys. Rev. Lett. 40, 589 (1978), D. Sherring- 
ton, Phys. Rev. B22, 5553 (1980), D. Sherrington, J. 
Phys. C14, L371 (1981). 

[14] W. G. Unruh, Phys. Rev. Lett. 46, 1351 (1981). 

[15] Artificial Black Holes, edited by M.Novello, M. Visser, 
and G. Volovick (World Scientific, Singapore, 2002). 



[16] Quantum Analogues: From Phase Transitions to Black 

Holes and Cosmology, edited by R. Schiitzhold and 

W.G. Unruh, Springer Lecture Notes in Physics Vol. 718 

(Springer, Berlin Heidelberg, 2007). 
[17] L. ,J. Garay, J. R. Anglin, ,J. I. Cirac and P. ZoUer, Phys. 

Rev. Lett. 85, 4643 (2000); L. J. Garay, J. R. Anglin, ,J. 

I. Cirac and P. Zollcr, Phys. Rev. A63, 023611 (2001); C. 

Barcelo, S. Liberati, and M. Visser, Class. Quant. Grav. 

18, 1137 (2001); P. O. Fedichev and U. R. Fisher, Phys. 

Rev. Lett. 91, 240407 (2003); V. A. De Lorenci, and R. 

Klippert, Braz. J. Phys. 34, 1367 (2004); H. Nakano, 

Y. Kurita, K. Ogawa, and C. Moon Yoo, Phys. Rev. 

D71, 084006 (2005); R. Schiitzhold, Phys. Rev. Lett. 

97, 190405 (2006); S. Wuster and C. M. Savage, Phys. 

Rev. A76, 013608 (2007). 
[18] T. Matsubara, Prog. Theor. Phys. 55, 351 (1955); R. 

Kubo, J. Phys. Soc. Jap. 12, 570 (1957); P. C. Martin 

and J. Schwinger, Phys. Rev. 115, 1342 (1959). 
[19] L. H. Ford and N. F. Svaiter, Phys. Rev. D51, 6981 

(1995). 

[20] K. Symanzik, Nucl. Phys. B190, 1 (1981); C. D. Fosco 
and N. F. Svaiter, J. Math. Phys. 42, 5185, (2001); R. B. 
Rodrigues and N. F. Svaiter, Physica A328, 466 (2003); 
R. B. Rodrigues and N. F. Svaiter, Physica A342, 529 
(2004); M. I. Caicedo and N. F. Svaiter, J. Math. Phys. 
45, 179 (2004); N. F. Svaiter, J. Math. Phys. 45, 4524 
(2004); M. Aparicio Alcalde, G. F. Hidalgo, and N. F. 
Svaiter, J. Math. Phys. 47, 052303 (2006). 

[21] E. Arias, N. F. Svaiter, and G. Menezes, Phys. Rev. D82, 
045001 (2010). 

[22] J. -P. Hansen and I. R. MacDonald, Theory of Simple Liq- 
uids (Academic Press, London, 2008). 

[23] P. W. Anderson, Phys. Rev. 109, 1492 (1958). 

[24] G. Krein, G. Menezes, T. C. de Aguiar, E. Arias and 
N. F. Svaiter Bose-Einstein Condensation in Inhomoge- 
neous Random Media, in preparation. 

[25] S. Fagnocchi, S. Finazzi, S. Liberati, M. Kormos, 
and A. Trombettoni, New J. Phys. 12, 095012 (2010) 
[arXiv: 1001. 1044 [gr-qc]]. 



